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Abstract 

r-| . Based on Regge-inspired arguments and counting rules, we formulate empirical models for 
zero-skewness generalized parton distributions (GPDs) H'^^^ and -E^^-* in the iso-vector sector. If 

Ph" 

D ] a hypothetical ai/p2 master trajectory is taken into account, we find that the polarized deep 
inelastic structure function gl , axial- vector form factor, pseudoscalar form factor, and lattice 
data are well described. Thereby, we use a symmetric valence scenario in which the 'spin puzzle' 
in the iso-vector sector is trivially resolved. Utilizing a minimalist 'holographic' GPD principle. 



> ■ 



lO ■ tying the t-channel angular momentum and collinear conformal spin together, we build skewness 
^ . 

O ' hard exclusive 7r+ electroproduction within the collinear factorization approach, we might conclude 

O ' that minimalist GPD models are disfavored at leading order. We provide then a set of GPDs on 



X 



dependent GPD models. Confronting these models with HERMES and JLAB measurements of 



the cross-over line that fairly describes both HERMES and JLAB measurements. 
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(a) (b) 

Figure 1: (a) Factorization of hard exclusive vr"*" electroproduction at large photon virtuality into a 
hard-scattering amplitude T, the flavor changing generalized parton distribution, and pion distribution 
amplitude (f)-,^. (b) Hard scattering amplitude at leading order. 

1 Introduction 

The hard exclusive meson electroproduction off proton, e.g., for vr^ 

e-{k)p{Pi)^e-{k')n{P2)'K^{q2), (1) 

is considered as a promising channel to study the partonic content of the nucleon. If the proton 
is probed with a longitudinal polarized photon with sufficiently large virtuality = —{k — k'Y, 
the amplitude factorizes in a hard part, calculable in perturbation theory, and non-perturbative 
quantities that are universal, i.e., process independent, however, conventionally defined pQ. As 
illustrated in Fig. [1] (a), the virtual photon knocks out a quark- antiquark pair from the nucleon. 
In the coUinear factorization framework the probability amplitude for this partonic process are the 
generalized parton distributions (GPDs), which depend on the longitudinal momentum fraction 
and a non-zero momentum flow in the t-channel [21 El S] • The formation of the meson is described 
within a distribution amplitude (DA) [5l[6], i.e., the minimal Fock state component of the meson 
wave function, where the transverse distance of the quark- antiquark pair is squeezed. 

The GPD framework is viewed as a powerful tool that provides insight in the partonic content of 
the nucleon and it ties perturbative and non-perturbative methods, for comprehensive reviews see 
Refs. [ZllH]. The hard exclusive processes, described to leading order (LO) accuracy, allow to access 
the GPDs on the cross-over line |9] . In particular, the hard exclusive photon electroproduction is 
widely considered as a theoretical clean process, see, e.g., Ref. [TU] and references therein. The hard 
exclusive meson electroproduction, which might be utilized as a flavor filter, was elaborated at LO 
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accuracy for numerous channels [TTl [T2| [T3| [H] and based on GPD models, estimates have been 
given, see, e.g., Refs. [151 [161 [17] . Perturbative next-to-leading order corrections were evaluated in 
Refs. [HI [ini [20] , and power suppressed contributions were studied in Refs. [211 [22l [23l [24] . Model 
studies show that perturbative corrections are rather large, see also the comprehensive studies in 
Ref. [25]. Moreover, the factorization breaks down for power suppressed contributions and so a 
systematic expansion does not exist [1] and various helicity amplitudes with leading twist-three 
suffer from so-called end-point singularities [21j . 

The primary goal of revealing GPDs from the hard exclusive meson electroproducion is rather 
challenging for both experiment and theory. First, the cross section parts, which arise from 
longitudinal and transverse polarized photon exchange, must be separated. To apply the coUinear 
factorization framework, the photon virtuality must be sufficiently large. Unfortunately, the onset 
of the perturbative framework is controversially discussed, where some considerations suggest that 
the perturbative description becomes reliable at a scale of > lOGeV^. Such a scale is larger 
than the virtuality that is accessible in present fixed target experiments. We like to emphasize 
that the emission and absorbtion of a quark from the target resembles the Feynman mechanism 
and intuitively this might be considered as well describable. On the other hand the formation 
of the meson can be affected to a larger extent by power suppressed contributions, related to its 
transverse size. Relying on this picture, improvements have been suggested that provide models 
rather than a systematic framework. 

For the hard exclusive electroproduction of a longitudinal polarized p\ meson it has been 
demonstrated that in a GPD framework, which contains also transverse degrees of freedom, ex- 
perimental measurements can be described [T5l [261 l27l [28] . In particular, at small xb the model 
[271 [28] , in which the GPD part is treated in the coUinear approximation, successfully predicts 
the spin density matrix. However, from the success of such models one cannot draw conclusions 
about the coUinear GPD approach, since model uncertainties are large and more importantly 
the interplay between perturbative and non-perturbative corrections in the coUinear factorization 
approach and the transverse degrees of freedom is essentially not understood. To judge on the 
coUinear factorization approach, one might confront theoretical "predictions" with experimental 
measurements. The problem of this approach is obvious, namely, one relies rather on a GPD 
model to reveal the GPD from experimental measurements. Thereby, the non-perturbative GPD 
models, not calculable from first principles, suffer from ad hoc assumptions, see, e.g., discussion 
in Ref. 

In this article we have a closer look to the phenomenological application of the GPD formal- 
ism to the hard exclusive vr^ electroproduction [TH [161 111]- More specifically, we consider the 
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longitudinal photoproduction 



72(gi)p(Pi) - n(P2)vr+(g2) (2) 

within the coUinear factorization approach at leading order of perturbation theory. We like to 
confront this approach with cross section measurements from the HERMES collaboration [30] 
and in HALL C at JLAB [HTl [32] . In the later experiment the longitudinal cross section could be 
separated, however, both the energy and the photon virtuality are rather low. At the HERMES ex- 
periment a Rosenbluth separation cannot be performed. Fortunately, according to Regge-inspired 
model calculations [33l[3ll[35], the pion pole enhances the longitudinal cross section at smaller val- 
ues of —t and it becomes the dominate part. We also include in our studies preliminary HERMES 
results for the single transverse proton spin asymmetry |36j . 

Our goal is to confront GPD models with measurements and to illuminate the interplay of 
model dependence and theoretical uncertainties. To set up our zero-skewness GPD models, we 
follow the idea, applied to polarized parton distribution functions (PDFs) to constrain or 
reveal non-perturbative quantities from generic arguments. Our framework is set up in Mellin 
space and it allows us to skew the resulting models in a convenient way, where we utilize the 
t-channel SO (3) partial wave (PW) expansion [38| [39l [7]. 

The outline of the paper is as follows. In Sect. [2] we recall the perturbative framework for 
hard exclusive vr"*" electroproduction at LO accuracy in the momentum fraction representation 
and then we present it in terms of conformal Mellin moments. In Sect. |3]we explore generic GPD 
modelling with Mellin moments, where we are guided by Regge-inspired arguments and counting 
rules. Thereby, we have to face the old problem that unnatural parity exchanges challenge Regge 
phenomenology. We will inspect the meson spectrum and we will argue that a odd J^~^ leading 
Regge trajectory might exist. As a side remark we shortly discuss the so-called 'spin puzzle' and 
we employ a symmetric valence matching scheme of collective and partonic degrees of freedom. 
Our zero-skewness GPD models are then introduced and we compare them with experimental 
measurements and lattice data. Relying on an ad hoc assumption, we set up a minimalist GPD 
model and inspect its properties. In Sect. H] we shortly recall the illness of the perturbative 
framework and we utilize a scale setting prescription that might cure it. We provide then model 
predictions for the longitudinal cross section and single transverse proton spin asymmetry and 
confront them with HERMES and JLAB measurements. We will then take the point of view that 
an improved LO formalism is reliable and we present GPDs on the cross-over line that describe 
experimental data. Finally, we summarize and give our conclusions. 
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2 Collinear factorization approach at LO 



We consider the hard exclusive 7r+ electroproduction off a transverse polarized proton in the 
target rest-frame. More specifically, the 2;-axis is directed along the momentum of the virtual 
photon and both the azimuthal angles (f) of the pion momentum and of the proton polarization 
vector are specified by the so-called Trento convention |30]. The cross section is straightforwardly 
evaluated in the collinear factorization approach to LO accuracy [T] . Thereby, the part of the cross 
section that arises from the exchange of a longitudinal polarized virtual photon is counted as the 
leading twist-two contribution. The exchange of a transverse polarized photon is on amplitude 
level formally suppressed by 1/ Q and will not be considered. 

In the following we can safely neglect the pion mass and mass difference of proton and neutron, 
i.e., we set M = Mp = M„. In the twist expansion of the cross section we will adopt here the 
same conventions as used for the hard exclusive photon electroproduction in Ref. ^1]. Namely, 
we decompose first the amplitude by inserting photon helicity states into a leptonic and hadronic 
part. We exactly calculate the leptonic part and drop any power suppressed contribution in 
the hadronic part. Thereby, we choose 1/Q as expansion parameter and the Bjorken variable 
xb = Q^/ (2gi - Pi) as scaling variable. Finally, we can write the differential cross section for hard 
exclusive vr^ electroproduction, induced by longitudinal polarized photon exchange, at twist-two 
level as 

da'^^ Tw-2 a^jjj sy'^ xb .2Mff 12 ^ 



dQ^dxBdtdip 1 - s ^/T~^^ 



-2 sin(0 - <Ps) ^1—^ e 53m ( H.+E:^_ 



"mm 



Here y = qi ■ Pi/k ■ Pi is the photon energy momentum fraction, t = {P2 — PiY is the squared 
momentum transfer, the symmetric scaling variable ^ is consequently replaced by xb/(2 — xb), 
and e = x-qM/ Q is a shorthand. The polarization parameter 



1 — y — hv'^e^ 



is the ratio of longitudinal and transverse photon flux. The minimal value of —t is 



"''^"^ 4xB(l-XB)+e2 ^1-xb' 

where the indicated approximation is valid for x^M"^/ <^ 1 — xb. The hadronic part is contained 
in the transition form factors (TFFs) and . They might be viewed as physical quantities 
that depend on the symmetric scaling variable ^, t, and Q^. 
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As depicted in Fig. [T^, the TFF (E^+) factorizes into a convolution of a hard amphtude 
Tud with the GPD H""^ {E""^) and the pion DA 

du J rfx 0^(M,/i^)r„d QV/^^ SV/^r) S }{x,C.,t\fi'^) + . . . , (6) 

where the elhpsis stands for power suppressed contributions. Here, we distinguish between the 
renormahzation scale fir and the factorization scale /i. Note that the truncation of the perturba- 
tion series induces residual scale dependence. The H^"^ and E'^'^ GPDs, depending on the quark 
momentum fraction x, ^, t, and /i, are flavor non-diagonal and parity-odd. They are defined by 
a form factor decomposition of light-ray operator matrix elements in the standard notation, see, 
e.g., Ref. [18]. The leading twist-two pion DA (p^^, depending on the quark momentum fraction u 
and factorization scale /i, is conventionally normalized to the pion decay constant: 

/ duMu,fi^) = U, /^^ 133 MeV. (7) 
Jo 

The short-distance dynamics, i.e., the hard amplitude Tud, can be perturbatively calculated 
as a series in the strong coupling as{fir), known to NLO accuracy [18]. From the tree diagrams, 
presented in Fig. [T)d, we recover the well known result [HI [161 fT5] : 



T.,Kx,e|«.(/i^), QVfi\ QVfil) = ^aM) \t. - ■n^ 1^ + ^(«' 



iVc i (1 - n) - a; - zO) n + a; - zO) 



si ' 



where Cp = 4/3, Nc = 3, and the fractional quark charges are Qu = 2/3 and Qd = —1/3. The 
hard amplitude ([8]) is proportional to the running couphng 

4:71 

c^si/jl) = — with /?o = 11 - 2nf/3 = 9 , (9) 

"qcd 

where we set Uf, the number of active quarks, to three and in the following we will take the QCD 
scaling parameter Aqcd = 200 MeV. We emphasize that in principle the renormahzation scale 
/iR can be ambiguously chosen because this ambiguity is annulled in the full perturbative series. 

In the following we will stay at LO accuracy and equate the factorization scale with the 
photon virtuality, i.e., we set /i^ = Q^. In the LO approximation the hard amplitude (IHl) possesses 
a factorized u and ^ =F 3; momentum fraction dependence!. Since the pion DA has definite charge 
parity-even, possessing the symmetry ^^^(u) = (f).^{l — u) , the convolution integral 

,2, X f\ Mu,Q'' 



UQ') = — (10) 



^Note that although NLO corrections might be rather large, it is expected that this factorization is only slightly 
spoiled by radiative NLO corrections, see, e.g., Ref. [H]. 
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enters as overall factor in the normalization of the TFFs ([6]). The value of the inverse moment 
ffTOj) is directly 'measurable' at LO accuracy in the fusion process 7*7 — > 7r° [5]. The factorized 
GPD part of the H"'^ and E'"^ TFFs ([6]) might be decomposed in terms with definite charge parity, 

■ ■ ■ ) = H^'^x, ■ ■ ■ ) ± ^"'^(-x, ■ ■ ■ ) , E'^'^^ix, ■■■)= E^'^x, ■ ■ ■ ) ± ^""^(-x, ■ ■ ■ ) • (11) 

Here the ± superscript refers in the first place to the symmetry behavior under x — —x and it 
coincides for both of our parity-odd GPDs with the t-channel charge parity C = ±1. Employing 
isospin symmetry, one can express the flavor off-diagonal GPDs to the common flavor diagonal 
ones in the iso-vector sector 



H^d ^ ^(3) ^Jjn_jjd^ ^ud ^ ^(3) ^ _ ^d _ ^^2) 

The convolution of these GPDs with the LO hard scattering kernel 

^(^,t, Q^) ^=5 f'dx ^^"^'^'^'^'^ , ^(F) = {n-''^{H^''^),r^^{E^^^)] (13) 
J ~i q — X — lU I J 

might be denoted as 'Compton form factors' ('CFFs'). Note, however, that only the charge even 
part (plus superscript) enters in the LO approximation of the deeply virtual Compton scattering 
(DVCS) amplitude. In the case of hard exclusive vr"*" production we have the combination 




^{3)+ + 1^(3)- 



(14) 




in which also the charge even part dominates. Putting all definitions together, we finally have a 
very compact LO expression for the TFFs 

(e,t, Q') =^ a.if^y-fuQ') I ~3)^ I ^~3), k,t, Q^) + • • ■ . (15) 

Here the renormalization scale is ambiguous and so <ys{f^r) niight be loosely viewed as an "effective 
charge". As we will argue below in Sect. HI its value might be extracted from the measurement of 
the pion form factor. Let us emphasize that in principle Eq. (fT5i) can serve as a qualitative test of 
the coUinear factorization approach, where the TFFs of hard exclusive electroproduction are 
expressed in terms of quantities that can be addressed in the LO approximation of the pion-to- 
photon transition form factor, pion form factor, and to a certain extent in DVCS. 

Let us also introduce the conformal PW expansion of the TFFs. The underlying symmetry is 
the collinear conformal group SL(2,M) and in pQCD it can be used to diagonalize the evolution 
equation, see, e.g., Ref. [13]. The conformal PW expansion for the pion DA is well known at LO 

EE], 



00 



V, Q') = V 6(1 - v)vCl/\2v - 1) E„,(Q, Qo)a^{Ql) , (16) 



even 
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where the conformal PWs are given in terms of orthogonal Gegenbauer polynomials of order 
n and index u = 3/2. The partonic "quantum number", conjugated to the momentum fraction u 
is the (integral) conformal spin n + 2. The first coefficient = 1 is fixed by the normalization ([7]) 
and does not evolve, i.e., Eo(Q, Qq) = 1. If the scale Q> Qo increases, the evolution operator 

reduces the magnitude of conformal moments a„ with n > 0. Here the harmonic sum is defined 
as Sm = YlT=i i- The inverse moment, given in Eq. (ITOl) as convolution integral, might be now 
expressed as a series of conformal moments 

oo 

lAQ') = 5^E„(Q, Qo)an{Ql) . (18) 

even 

For the t-channel crossed GPDs, i.e., the so-called generalized DAs, one might write down an 
analogous conformal PW expansion as in Eq. (fT6!) . The conformal moments of a quark GPD are 
group theoretically defined in terms of Gegenbauer polynomials Cj of order j and with index 
u = 3/2: 

iv,t, Q') = 2!!fr|3)2 + ^] //-^^c^f (-/^) |f } (-,^,^, Q') ■ (19) 

Note that our definition ensures that [anti] symmetric functions possesses only even [odd] integral 
moments, which are even polynomials in rj of order j [or J — 1]. In the limit ?7 the conformal 
moments of H GPD coincide with the Mellin moments of a polarized PDF: 

Hf{r^ = 0,t, Q2) = j\xx^Hf{x,r] = 0,t, Q^) . (20) 

As in the case of a SO (3) PW expansion, labelled by the t-channel angular momentum, such a 
series does not converge for the s-channel analog. We can borrow us the known recipes and employ 
a Sommerfeld- Watson transform [H], which reads for the 'CFFs' as follows 

(f llo 1 p,, J|^,,/2)-.-.E,(Q.Q.)il^(})K,,SS). (21) 
[^J 2tJ^_i^ r(3/2)r(3 + j) smvrj ]^Ejj 

The analytic continuation of the evolution operator is straightforward by replacing n ^ j and 
representing the harmonic sum Sj+i = + 2)—iIj{1) in terms of Euler's ip functions. The analytic 
continuation of integral conformal GPD Mellin moments (fT9l) . is dictated by the requirement that 
after crossing ,^ ^ 1/,^, the Mellin-Barnes integral yields a series, rather similar to the convolution 
integral (ITS]) . In general this requires a separate continuation of even and odd conformal moments. 
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3 Generic modelling in Mellin space 



A GPD F{x, rj, t) at a given input scale /i^ is an intricate function of three variables, the s- 
channel momentum fraction x, its t-channel counterpart rj ~ (Pi — ^2)+; and t. With the growing 
amount of experimental data for the hard exclusive photon electroproduction it became obvious 
that popular ad hoc GPD models do not posses predictive power on quantitative level; a more 
detailed discussion can be found in Ref . [21] . For our purpose, namely, a qualitative confrontation 
of the coUinear factorization approach with hard exclusive tt^ electroproduction measurements, it 
is sufficient to deal with zero-skewness GPD models and an ad hoc prescription to skew them. In 
the following we work with conformal GPD Mellin moments (fT9l) . which in principle allows us to 
set up flexible GPD models. 

To parameterize the ?7-dependence of GPDs, we follow Ref. [SB] and expand the conformal 
moments (fT9l) in terms of t-channel SO (3) PWs. These PWs are given in terms of the Wigner 
rotation matrices (ioo(cos6') and dQi{cos6), i.e., by Legendre or Gegenbauer polynomials, re- 
spectively [U US]. Here the second subscript refers to the helicity difference of the final state 
jSfUDpjl in the t-channel center-of-mass frame, i.e., for (anti-) aligned nucleon helicities it is zero 
(one). The center-of-mass scattering angle 6 is approximatively related to the skewness parameter 
7] ^ —1/ COS 6. Note that the exact relation would yield more cumbersome formulae that provide 
no advantages or improvements for our GPD modelling. The crossed version of the Wigner rota- 
tion matrices are noted for shortness as d-QQ^f]) and c^qiI'?) and they are normalized in the limit 
?7 ^ to one |15|. Their analytic continuation with respect to the t-channel angular momentum 
J is defined in the traditional way and we express them by hypergeometric functions: 

?7 — 1"^ 

2-^1(1/2 + J) ■' 2 

2-^1(1/2 + J) ■' 1 
Now we can straightforwardly write down the S0(3)-PW expansion for integral conformal 
GPD moments (fTOl) . For the target helicity conserved GPD H the S0(3)-PW expansion reads 

HtM = E ^ ^ ^~^^' Hfj_,{t)rj^+'~'diM . (24) 
j=i 

Here even j moments (plus superscript) and odd j moments (minus superscript) are build up by 
odd J~^^ and even J exchanges, respectively [391 17]: 

= 3++, ■ ■ ■ , (j + 1)++ , = 2-, 4--, ■ ■ ■ , (j + 1)- . (25) 

Note that here the difference j — J is always odd and so a confusion with the notion signature, which 
might be used in both the conformal and SO (3) PW expansion, could appear. We consequently 



^fM-'4m^^^^F^(-'-'.^''^)- (23) 
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stick here and in the following to the Regge terminology and associate signature plus and minus 
with even and odd J, respectively. 

The t-channel helicity conserved GPD moments are expressed by a combination of s-channel 
flip and non-flip contributions: 



' i±r-i)^ 



J=0 



t 



ri'-' 'iUri) . (26) 



Now the even j moments (plus superscript) and odd j moments (minus superscript) are associated 
with even J ^ and odd J+~ t-channel exchanges, respectively [391 17]: 

J^^ = o"^2"^•■■,r^ j^^ = i+^,3+-,---,j+-. (2?) 

We might utilize the crossed version of Wigner's rotation matrices (ig ol'?) within the assignment 
(|27|) in the SO (3) PW expansion of E GPD conformal moments: 

E}{V, t) = Y. ^-^^^EfA^) V'-' dioiv) ■ (28) 
j=o 

Then the H GPD conformal moments can be alternatively represented by the SO (3) PW expansion 
and fl26|l in terms of the amplitudes Hjj^i(t) and H'-j{t), respectively. Consequently, these 
both PW amplitudes can be mapped to each other. Below we will assume that the Hj conformal 
moments are only built by one leading S0(3) PW J = j + 1 with odd or even J . This 
implies that all PW amplitudes of Hjj(t), appearing in Eq. (l26l) and containing contributions from 
even J""*" or odd J~^~ exchanges, are exited. Of course, they can be entirely expressed in terms of 
Hjj(t) amplitudes: 

(t) = r(j + 3/2)r(j + i) ~ j>j>o, (29) 

jjK J 2-^-:'r(j + 3/2)r(J+ 1) nKJ^ J- - V ; 

In the following sections we model the zero-skewness Mellin moments, 

Hfiv = 0,t) = H^fit) and Ef\v = 0,t) = ^]?(t), 

in the iso-vector sector. Although Regge theory is developed for on-shell scattering, it is obvious 
from the phenomenology that effective Regge behavior also appears in off-shell amplitudes. As a 
fact of matter, in the coUinear factorization approach effective Regge behavior is encoded in PDFs, 
obtained from global fits, and so also in GPD models. Moreover, it is rather robust under evolution, 
except for the "pomeron" -like exchange [16]. Hence, Regge poles should be a building block for our 
partonic SO (3) PW amplitudes. The corresponding PW amplitudes will be interpreted according 
to Regge phenomenology and they are in the t-channel, i.e., 72 tt" — ^ np, given by exchanges with 
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definite J. Analogous as in Regge theory, applied to high-energy on-shell scattering, this allows us 
to resum t-channel exchanges and, finally, they determine the small-x behavior of GPDs or PDFs. 
Moreover, a link between effective Regge behavior and t-dependence of nucleon form factors has 
been also employed or observed in zero-skewness GPD models [171 HHl SI]. We are going along the 
line of Ref. [HJ HH HQ] and unpublished work, shortly mentioned in Ref. [50] (see corresponding 
talk at DIS2008). There the t-dependence of H and E GPD Mellin moments is modelled as a 
product of leading and daughter Regge poles. Since the empirical H and E GPD models describe 
within generic numbers, obtained from counting rules and Regge theory, form factors, unpolarized 
PDFs, and lattice data, we are encouraged to utilize generic modelling in the following sections. 

However, first we have a look to the phenomenological status of Regge trajectories in the 
meson spectrum. We then build our zero-skewness GPD models and confront them with the 
polarized deep inelastic scattering (DIS) structure function gi, form factor measurements, and 
lattice data. Finally, we skew our GPD models within a minimalist prescription as they arise 
from the dominance of the leading SO (3) PW in the expansion (124]) and (1281) . We also discuss the 
model properties, relevant for the phenomenology of hard exclusive vr"*" electroproduction. 

3.1 Regge trajectories from meson spectroscopy 

According to the quantum numbers fl2^ and (1271) . we are dealing with unnatural parity 

p = -{-iy 

t-channel exchanges such as TfQ^/bf^ and 0'i'^/p2~ ones. Such exchanges challenge Regge phe- 
nomenology, applied to high energy on-shell scattering processes, see reviews [S21IS3]- We are not 
aware that the problems there, e.g., an unexpected large ai pole contribution in 7r~p p^n [52] . 
were finally solved, e.g., in terms of very large Regge cut contributions. Therefore, and because 
of the immense experimental progress in meson spectroscopy during the last decade, see the com- 
prehensive reviews [M| it is worth to have a fresh look to Regge trajectories in the meson 
spectrum. Thereby, we take experimental data from the Particle Data Group (PDG) [51], where 
we also include mesons listed under further states. 

According to Eq. the leading SO (3) PW amplitudes H^-f^ and H!jf~ with even and odd 
j = J — 1, respectively, carry odd J^~^ and even J t-channel quantum numbers. We associate 
them with meson exchanges that belong to ai and p2 trajectories. In the Chew-Frautschi plot in 
Fig. [2] the iso- vector mesons 

ai, as, p2 and p^ 

might lie on trajectories, which we traditionally consider as linear equidistantly spaced, i.e., 

J(M2) = (M^) , a W (M^) - ao - ^ + a'M^ , n = 0,l,2,---, (30) 
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^ p/tt 


0.50 + 0.87 M^/GeV^ 




aalp = 


0.53-n + O.S^M^/GeV^ 




a„ib= 


-0.014 + 0.72M^IG&V^ 






nonlinear 




M^(GeV^) 



Figure 2: Spectrum of iso- vector mesons oi, 03, p2 and /04, the states of the leading vro/61 and pi/a2 
trajectories. Curves show linear Regge trajectory fits, ai/p2 (dashed), vro/fti (dot-dashed) and pi/a2 
(solid), and an example of a non-linear vro/fti trajectory (dotted). Spectroscopic data are from the 
Review of Particle Physics [51] . 



where a(M^ 



a 



(n.=0) 



(M^) is called master trajectory. In particular, the states 



ai(1260)[1230±40], 03(2070) [2070 ± 30], p4(2240)[2230 ± 25] 

[values in square brackets are the mass and its error in MeV, quoted by the PDG] might belong 
to the trajectory with the largest intercept, called Oi trajectory. Note however, that p2 is not 
established and 

a3(1875)[1874±43±96] 

lies between the pi/a2 master and the ai trajectory within an equal distance of slightly more than 
1-cr standard deviation (errors are added in quadrature). The states ai(1640), p2(1940), a3(2310) 
and ai(2095), p2(2240) might belong to two daughter trajectories. 

Let us first suppose that both states a3(1875) and a3(2070) nearly belong to the ai trajectory. A 
least square fit within a linear ai trajectory and two daughters reveals that x^/d.o.f. ~ 4.1/(9 — 2) 
does not contradict our hypothesis and the result 

-(0.32 ± 0.30) + (0.83 ± 0.07) s/GeV^ , xVd-o.f. ^ 0.59 (31) 
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confirms the ai trajectory that is "established" in the hterature. The trajectory fl3T]) has a typical 
t-slope and its intercept —0.6 < ttai(O) < suffers from a large uncertainty. Such an uncertainty 
appears also in the analyzes of high-energy scattering processes, see, e.g., discussion in Ref. |56j . 

Certainly, other interpretations of the aj meson spectrum are not excluded, e.g., one might 
consider the Oi trajectory as degenerated with a non-linear ttq trajectory [5^, as illustrated by 
the dotted curve in Fig. [2l Also it has been argued that chiral symmetry might be partially 
restored at large J, and so exotic parity doubling^ appears in the meson spectrum, see review [58] . 
According to that, at sufficiently large J a meson with natural parity, belonging to the master 
trajectory, should posses two partners with unnatural parity. The author of Ref. [59] argued that 
the P3~(1690) meson has a 3"''^(~ 1700) partner. The large mass error allows us to speculate that 
a^'^i^lSTd), reported in 2002 in Ref. [60], is a parity doubling candidate, which might indicate a 
partial chiral symmetry restoration. 

We now count the a^~^{1875) state to a master trajectory and like to know whether it is 
degenerated with the pi/a2 trajectory. A fit to the ai/p2 nieson spectrum yields the result 

= (0.53±0.34) -72 + (0.86±0.08) s/GeV^ xVd-o.f . ^ 0.67 , (32) 

which is with n = 1 fully compatible with the previous one (13T|) . Although the uncertainty is large, 
the means somehow suggest that this trajectory is degenerated with the pi/a2 one. In Fig. [2l our 
fit result ( 1321) . shown as dashed curve, is hardly to distinguish from the pi/a2 master trajectory 
(solid curve), which we obtained from a separate fit to its first five members. 

Further support for the existence of a master ai/p2 trajectory arises from polarized DIS. Here, 
the rise of the iso- vector polarized nucleon structure function g^\xB) at small-xe might be viewed 
as a Xg'^ power behavior with A ~ 1/2 [611 ESj, illustrated below in the left panel of Fig. [3l Such 
a power behavior, partly even a stronger one, is included in standard PDF parameterizations and 
it is predicted by a double log resummation [63] . On the other hand in the Regge framework, one 
equates A = aai(O) [64] • Hence, from the "established" ai trajectory (!3T]) it is then expected that 
9i^\xb) in the small-xe region vanishes or approaches a constant, which contradicts experimental 
measurements. As the authors of Refs. W2\ , we are not aware that this Regge puzzle is solved. 

Although it is not conclusive from the meson spectrum that a leading trajectory for unnatural 
parity states exists, we conjecture a master ai/p2-trajectory, 

aa,/p^ (s) ^ ap,/a^ (s) = (0.5 ± 0.05) + (0.87 ± 0.05) s/GeV^ , xV^-o.f . = 0.24 , (33) 



^ Using J = q;'^")(M^) and the assignment P = (—1)^+^ and C = (—1)^+'^, where according to the quark model 
S — {0, 1}, it is easy to see that a parity doubhng pattern appears in the meson spectrum, except for the master 
trajectory ap^^g^^{M'^) itself. Partial restoration of chiral symmetry predicts that parity (— 1)"^ states on the pi/a2 
master trajectory at sufhciently large J are degenerated with two further parity (—1)"'+^ states, carrying C = +1 
and C = — 1; respectively. 
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to explain the rise of the polarized DIS structure function gl . We consider this hypothetic Regge 
trajectory as degenerated with the master pi/a2-trajectory and we view the "established" ai- 
trajectory as its first daughter. We did not investigate whether our conjecture resolves problems 
or contradicts phenomena in on-shell processes at high energy [52l [53] . 

We inspect now the even J""^ and odd J~^~ Regge trajectories that appear in the SO (3) PWs 
of E^j and E~j, respectively, where j = J, see Eq. ( 1271) . Presumably, the leading trajectories 
are non-linear and their daughters are populated with natural parity vector states vrf^, having 
exotic quark model quantum numbers (5 = and even L). Since this is not compatible with 
our hypothesis fl5Ul) . we do not utilize daughter trajectories. The leading signature even and odd 
trajectory are populated with 

71^ [(139.57018 ± 0.00035) MeV] , 7r2(1670) [(1672.4 ± 3.2) MeV] , 7r4(2250) [(2250 ± 15) MeV] 
and 

6i(1235) [(1229.5 ± 3.2) MeV] , 63(2025) [(2032 ± 12) MeV] 
mesons, respectively. 

If we ignore the 7r4(2250) state, reported in 2001 in Ref. [53], and we assume that the trajectories 
are linear and degenerated, we find from a fit the 7ro/6i-trajectory 

ano/bAs) = (-0.014 ± 0.002) + (0.72 ± 0.01) s/GeV^ , xV^-o.f. ^ 1.1 , (34) 

shown in Fig.[2]as dot-dashed line. Since this trajectory is often used, we denote it as "established". 
Compared to a typical meson Regge trajectory, the "established" vro/fci-trajectory has a smaller 
slope parameter and much smaller intercept, which is close to zero. Remarkably, the intercept can 
be numerically expressed as a7r(0) = — m^a^. This allows us to parameterize the pion trajectory 
in terms of the pion mass and the slope parameter: 

(s) = -a' (ml - s) with = ^ 0.14 GeV. (35) 

V < 

However, the "established" trajectory (IMl) is disfavored by the fifth state 7r4(2250) (a fit within 
a linear trajectory to all five states provides x^/d-o.f. ~ 3.5). From Fig. [2] it seems to be obvious 
that the anomalous behavior of the trajectory appears at low J values rather than large ones, see 
for example curve (dotted) and compare also with master trajectory. If we now adapt the heavy 
pion world (m^ ~ 600 MeV^), the vro/61 trajectory (l35l) could be viewed as linear and parallel to 
the master one. On the other hand in the chiral limit — s> the pion might be viewed as a 
singular state that does not belong to a linear 6i/7r2-trajectory, connecting the other four states. 

Presumably, the non-linearity of the vro/foi-trajectory is related to spontaneous chiral symmetry 
breaking. However, we believe that the chiral limit does not provide a perfect guidance to model 
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builders that are interested to describe the real world. Thus, we should in fact utilize a non-linear 
pion trajectory rather than a pion pole. However, for simplicity we compromise and choose the 
hypothetic a\l p^- and vro/^i-trajectory to be linear within a common slope parameter that is 
smaller and larger than those of the "established" ai/p2- and vTo/fci-trajectory, respectively, 

1 4 

a{t) = a + a't , a = - , a' = g ! cunit) = a.„ + at , = —m^a , = 0.14 GeV . (36) 

2 5 GeV 

Here and in the following we denote a leading master trajectory as a{t) and its intercept with the 
shorthand 

a = a{t = 0). 

3.2 Modelling H^f 

3.2.1 t = case: polarized PDFs 

We consider first the t = case of the leading SO (3) PW amplitudes in conformal GPD moments. 
They coincide then with the common Mellin moments of the polarized PDF: 

fijf«^o).Ar for .-{:;;|:;. 

cf. Eq. ( 120|) . For even and odd j the parity charge-even and -odd type PDFs are to be taken 
Aqf^ = dx Ag(3)^ (x) , Ag(3)± = Au - Ad± {Au - Ad) , j = | ^ J^"" . (38) 

We also distinguish between valence, sea, and anti-quarks, where the sea content of polarized 
quarks is expressed by the anti-quarks Aq{x). The polarized charge parity-even PDF is then 
given as sum of valence quarks and antiquarks and the charge parity-odd one is equated to valence 
quarkj^: 

Ag+(x) = Ag™'(|x|) + Ag^°"(|x|), Ag^^"(|x|) = 2Ag(|x|) , Ag-(x) = sign(x) Ag™i(|x|) . (39) 

Let us make a side remark with respect to the 'spin puzzle', which we understand as the mis- 
match of quark model and partonic degrees of freedom. In the quark model the nucleon possesses 
no antiquarks and its constituents are considered as collective degrees of freedom. Loosely spoken, 
the constituent quarks are somehow made up from the partonic content. If one literally identifies 
the partonic and quark model degrees of freedom, one easily generates momentum fraction or spin 
'puzzles'. It might be somehow legitimate to identify the collective constituent quarks with the 
partonic valence quarks [HS] and one might wonder whether the 'spin puzzle' exists then. The 'spin 

^ For convenience, we proceed in momentum fraction space and below in Sect. 13.2.21 we return to Mellin space. 
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puzzle' is perhaps easily solved by an appropriate matching prescription of collective and partonic 
degrees of freedom or, as in the case of the 'momentum fraction puzzle', one simply accept that 
one deals with different degrees of freedom. Of course, the decomposition of the proton spin in 
terms of its fundamental degrees of freedom remains a serious challenge. 

Although it is rather likely that SU(6) symmetry in polarized valence quark PDFs holds true 
at small-x values and it is broken in the large-x region, we find it interesting to explore the valence 
matching scheme of collective and partonic degrees of freedom. To do so, we adapt flavor-spin 
SU(6) [SU(3) ® SU(2)] symmetry to fix the normalization of polarized valence quarks in the 
iso- vector sector by the group theoretical value of the axial- vector coupling g^^: 



Jo 



val 



suj) 4 _ -1 ^ 5 

3 3 3 ^ ^ 

The axial- vector charge qa, i.e., the expectation value of the axial- vector current in the forward 
kinematics, is given by the lowest moment of the polarized charge parity-even quark PDFs. Hence, 
both valence and sea quarks contribute: 



9A 



/ Ag(3)+(a;) = M^'™' - M^'™' - rfT'™' + c/^'™' + 2{u^ -u^ -d^ + d^' 

'^i' 1.2670 ±0.0035, (41) 



and the experimental value, taken from Ref. [Slj, is smaller than the SU(6) quark model expecta- 
tion. For an isospin symmetric polarized sea, the axial- vector charge qa reduces to its valence part 
5f™^ In such an ad hoc scenario, historically used for simplicity, e.g., in standard PDF fits, the 
SU(6) value contradicts the experimental finding fHT]) . We clearly spell out that the SU(6) 
value fHUl) is not ruled out and a 'spin puzzle' does not exist. Rather SU(6) symmetry forces us 
to introduce an isospin asymmetric sea with a large negative net contribution 

^^^'^ = j Ag(=^)^<=^(x) = 2{u^ -u^ -d^ ^'d}^ 

suj)_2_ ^^2) 

Here we used that qa = 5/3 — 2/5 ~ 1.267 is the mean of the experimental value (14T1) . As 
we will see below, this SU(6) symmetric valence scenario is also not ruled out by polarized DIS 
measurements and some lattice simulations. Our scenario is in conflict with chiral quark soliton 
model (xQSM) considerations [671 [68l [69l [70| [7T| . predicting a flavor asymmetric sea within a 
positive value for g^^"^] however, qualitatively supported by meson cloud model considerations, 
suggesting a (small) negative value of g^^^ [721 ESI EH 175] . 
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The next ingredient for our generic PDF model are large-x counting rules, obtained from 
dimensional counting or diagrammatic analysis [76| [77]. They state that for a proton with helicity 
+1/2 a quark PDF vanishes for x — ^ 1 as (1 — x)'^, where the power 

7 = 2n, - 1 + (1 - 2h) (43) 

is determined by the number of spectators Ug and the helicity h = ±1/2 of the struck quark [37] . 
For valence and sea quarks with helicity aligned to the proton one the exponent 7 is given by 
(3 = 3 and (3 = 7, respectively. If the struck quark helicity is anti-aligned these values increase 
hj 6(3 = 2 to (3 + 6(3. We recall that the small-x behavior of our polarized PDF is determined by 
the intercept of the hypothetical master trajectory fl33p . i.e., \imrc^QAq{x) oc which we set 
to the generic value a = 1/2 . 

The small and large-x behavior are incorporated in a factorized ansatz: 

Ag(3) (a;) = ^(3) (a;^ ^ = 0, t = 0) (44) 

(1 -x)^-/J(l -x)^+'^^ 
B{l-a,l + (3)-hB{l-a,l + (3 + 6(3) 

where the normalization of the lowest Mellin moment is given by the corresponding partonic axial- 
vector couplings (7™^ and g^^^. In the partonic model interpretation the h parameter adjusts the 
ratio of quarks with aligned and anti-aligned helicities. We might assume that at small-x the 
sea quark density becomes helicity independent, i.e., we set h^°^ = 1. For this specific value the 
polarized PDF 0441) behaves with a = 1/2 in the small-x region as = x^/^, which is nothing 
but the Regge behavior that belongs to the established ai trajectory, cf. Eq. (!3Til . [78] . However, 
the partonic interpretation of h^°^ might be too rigid for /i™^ Namely, at large-x also a certain 
amount of valence quarks with aligned helicities might be suppressed by some additional powers 
of (1 —x)^^ [37j. We choose A/5 = 6(3 = 2 and consider now /i™^ as an effective parameter. In the 
following considerations we set /i™' 0.85. Relying on our partonic inspiration, the hypothetical 
Regge intercept, SU(6) symmetric valence scenario, and canonical counting rules, we can fix the 
model parameters: 

a^=i, a'^=^, ^?r''^=i'^^, /i™'P=0.85, /?™' ^= 3 , 6^^' (3 ^=2, (45) 

^_su(6)_2^ /j-as-i^ /^--l^'^T, 5^'^^/?^=2, 
5 

where also the generic Regge slope parameter, see Eq. (136|) . is included. 

A few comments are in order. For our convenience we take the parameters, specified in (l45l) . 
at the input scale /i^ = 2GeV^. The phenomenological values for (3 and 6(3 should in accordance 
with evolution be moderately larger while the hypothetical Regge intercept a is rather robust 



17 




Figure 3: Standard LO PDF parameterizations, described in the text, and a generic PDF model (solid), 
specified in Eqs. (glD and (gS]), are confronted with HERMES results [79l [80] for the DIS structure 
function gl (x) (left) and (Au — Ad){x) (right), where only statistical errors are shown. 

under evolution. Instead of a SU(6) symmetric valence scenario we could also have adopted a 
standard one from the literature. After an appropriate change of /i™i/sca^ g^jj these scenarios 
allow us to describe the measurements, discussed in this article. We add that the model can be 
easily extended to unpolarized PDFs. In such a framework the parameters /i™'/sca^ considered 
as effective ones, can be constrained from the large-x behavior of DIS structure functions. It is 
beyond the scope of this paper to present such studies. 

The iso-vector part of the polarized structure function, expressed at LO by the polarized PDF, 
can be experimentally obtained from DIS measurements on proton and neutron target: 

9?\xB, Q') = g^^\xs, Q') - gt\xB. Q') (46) 
"i'W=^)+(x = XB,/i^ = Q^). 

Here we conventionally set the scale /i^ = Q^. Our generic PDF model fairly agrees with standard 
polarized PDF parameterizations, taken at the input scale jj? = 2 GeV^. This is illustrated in 
the left panel of Fig. [3l where our model (solid) is displayed together with measurements of 
the HERMES collaboration [79], and standard parameterizations: Gehrmann/Stirling A (95) [81] 
(dashed), Bliimlein/Bottcher (02) scenario 1 [82] (dotted), and Gliick/Reya/Stratmann/ Vogelsang 
(00) 'broken valence' scenario [83] (dot-dashed). Note that we fixed in all PDF parameterizations 
= 2 GeV and that the various parameterizations have a rather different partonic interpretation. 
The GS and BB parameterizations rely on a flavor symmetric polarized sea, the 'broken valence' 
GRSV scenario possesses a positive sign for the polarized sea in the iso-vector sector, contrarily, 
to our SU(6) symmetric valence scenario which requires a negative sign. Hence, the valence and 
sea quark content of the axial-vector charge is rather different, too. Also the small-x behavior 
of the various parameterizations differ. Since the areas under the PDF curves in the left panel 
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of Fig. [3] are fixed by the Bjorken sum rule (or equivalently the normahzation condition pTl) for 
the axial- vector current), the small-x behavior of the various parameterizations is essentially a 
consequence of the assumed large-x behavior and the specific parameterization, in particular, the 
flavor scheme for the polarized sea. 

Semi-inclusive pion production allows addressing polarized antiquarks in the iso- vector sector. 
In the right panel of Fig. [3] we show the partonic interpretation of a HERMES measurement 
|80j . obtained from a LO fit (full circles), together with various PDF parameterizations. The 
'measurement' suggests that x{Au — Ad){x) is highly non-symmetric, where its sign might even 
change in different x regions; however, statistically, the asymmetry is also compatibly with zero. 
The figure also illustrates that present experimental data do not allow to quantify the polarized 
sea and that even a qualitative judgement on the different scenarios can not be drawn. This 
conclusion is based on more solid ground, namely, on the extensive analysis of polarized PDF 
uncertainties in global fits from Ref. [H]. Our SU(6) symmetric valence scenario is supported by 
global fits from Ref. [85] . 

Because of the polarized PDF uncertainties, one might wonder whether the so-called 'spin 
puzzle' or 'crisis' in the iso-vector sector has ever existed. The variety in the behavior of standard 
polarized PDF parameterizations might be also rather important for GPD model builders. Namely, 
these PDF uncertainties might propagate into GPD model predictions for hard exclusive processes. 

3.2.2 t-decoration via Regge trajectories 

To decorate the conformal GPD moments (!37|) with t-dependence, we rely on dimensional counting 
rules for the large —t behavior of form factors. Lattice results clearly show that with increasing 
j, i.e., in the momentum fraction representation in the limit x —>■ 1, the t-dependence dies out 
[86| [87] . To unify these features with Regge behavior, we write the t-dependence as product of 
monopoles, build up from the leading Regge trajectory and its daughters: 



where the Pochhammer symbol for integral p > is defined as 

{x)p = x{x + 1) ■ ■ ■ {x + p — 1) 

and the intercept a has the same value as in the polarized PDF ( l44l) . Hence, the zeros of the 
pre-factor in our model (H7|) cancel the corresponding poles in the PDF moments, which are then 
replaced by a product of t-dependent Regge poles 1/[1 + j — (y(t)] ■ ■ - [p + j — a(t)]. The number 
p determines the asymptotic of the PW amplitudes ( l47l) at —t oo, which behave as (— t)"^. 
Employing the Drell-Yan exclusive-inclusive relation 7 = 2p — 1 [88] , we should choose p = 2 and 
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Figure 4: The left panel shows the axial- vector form factor versus —t: world fit ()48|) [error band], generic 
H GPD model ()47p [solid], as specified in the text, and lattice results (full circles). In the right panel 
the same GPD model is confronted with lattice data for generalized form factors: j = 1 (dot-dashed, 
triangles) and j = 2 (dashed, squares). The scale is = 4GeV^ and lattice measurements are taken 
from data set VI of Ref . |89] . 



p = 4 for valence (7 = 3) and sea quarks (7 = 7), respectively. To simplify the model procedure, 
we take for both contributions p = 2 and we complete the Regge intercept to a generic leading 
meson trajectory (H5l) . 

We illustrate in the left panel of Fig. H] that the t-dependence of the monopole product ansatz 
fjTr|) for j = 0, displayed as solid curve, is in fair agreement with world dipole fits (error band), 

gA{t) = — ^^^^^ 2 , = (1.069 ±0.016) GeV^ [(1.026 ±0.021) GeV^] , (48) 

to the axial-vector form factor measurements from tt"*" electroproduction and neutrino scattering 
experiments [M\ given in square brackets], taken from Ref. [5U]. We also display lattice data (full 
circles) from the set VI of Ref. [SH] , which at larger —t considerables overshoot the experimental 
measurements (only statistical errors are shown) and possess at low —t a rather flat t-dependence. 
The dotted curve is included for illustration and it arises from a shift of our solid model curve. In 
the right panel of Fig.lHwe plot the second (dot-dashed) and third (dashed) moment 0471) versus —t 
for fx^ = A GeV^, where the polarized PDF moments are specified in Eq. (1451) at /i^ = 2 GeV^. Our 
generic model predictions are now compatible with lattice data. We emphasize that in particular 
the lattice result for the j = 1 valence generalized form factor at t = 0, given aqj 0.25 ± 0.02, 
does not "rule out" the SU(6) symmetric valence scenario prediction ^ 0.29. In standard PDF 



''Lattice measurements suffer from systematical uncertainties that can hardly be estimated and so it is common 
to quote only the statistical error. 
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parameterizations the j = 1 Mellin moment is typically smaller, namely, ~ 0.2. Note that this 
discrepancy in the valence sector essentially originates from the difference between non-standard 
flavor scheme (l45l) and flavor symmetric polarized sea. 

Finally, we summarize our generic model. Analogous as for PDFs in Eqs. f l55]) and fl5Ul) . the 
odd j (or — ) moments are entirely associated with valence quarks while the even (or +) moments 
are given by the sum of sea and valence quark contributions. We plug the Mellin moments (138|) 
of our PDF ansatz (H4l) into Eq. (1471) and obtain so the leading SO (3) PW amplitudes 

^3)+^^3)val^^3)sea^ ^3)- ^ ^3)val ^ ^.^^ ^^g^ 

H(^\t)= B{l+j-a,l + (3)-hB{l+j-a,l + f3 + 6f3) 

" B{l-a,p + l)-hB{l-a,l + p + 6p) 

where the model parameters are specified in Eq. (H5l) . 



3.3 Modelling E^f 

Analogous as in the previous section, we build now a generic model for the leading SO (3) PW 
amplitude Ejj=j{t). Thereby, we borrow us the j-dependence at t = from the ansatz (HOj) for 
Hjj(t = 0) with generic Regge intercept a = 1/2. The new {/?, SP,e} parameter set {h is replaced 
by e, and to lighten the notation, we use the same symbols /5 and 6j3 as for H) will be specified 
below. 

At large —t the helicity non-conserved moments Ejj{t) should decrease by one power of (— 
faster than Hjj(t), i.e., it should generically behave as l/(— t)^. In the heavy pion world, in 
which the vro/6i trajectory is presumably degenerated with the ai/p2 one, we would write the 
i-dependence of Ejj(t) as a product of three monopoles. 



(where j = J) 



(j - a At)) 3 

which are expressed in terms of the iTo/bi Regge trajectory and its first two daughters. In the real 
world we choose the Regge trajectory (1351) . however, we assume that its two daughters are now 
degenerated with those of the generic master trajectory a{t). Finally, we introduce the partonic 
decomposition of even and odd Mellin moments and we write our ansatz for them as follows: 

^3)+ ^ ^3),val ^ ^3),sea ^ ^3)- ^ ^3),val ^ ^^q^ 

~(3) ^ <GeV^ (l+J-a)2 ^(1 +3-aA + (3)- eBjl + j - g^l + f] + Sf3) 
^ ^ ^^j-a^(t)(l+j-a(t))2 B{l-a,l + p)-eB{l-a,l + P + Sp) 
where the generic and pion Regge trajectories are specified in Eq. fl5B]) . respectively. Our model 
for even j is normalized at j = J = 0, 

E,, {t) = gp{t)-—^-^^—-^^, 9p = 9p-9p+9p, (51) 
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by the pseudoscalar coupling gp. In our partonic interpretation it contains a valence and sea quark 
part, where the former one also determines the normalization of the odd j (or — ) GPD moments. 
The lowest moment ( ISTl) is nothing but the pseudoscalar nucleon form factor. As designed, our 
parameterization fl51l) possesses the pion pole, which appears by means of Eq. fl5^ as follows: 

K j=0 Oi'-K 1 .,1, 2 "tt 

^ — With — 



j — a-„{t) — — a'^t — t ^ a' 

Let us show that our generic GPD model satisfy the well-known low energy constraints that 
arise from the hypothesis of a partial conserved axial- vector current (PCAC). PCAC relates the 
axial- vector and pseudo-scalar form factor, in our notation the GPD combination fl26|l for j = 0, 
with the pion-nucleon form factor: 

~«S.it)^i^M^)^^^^G,Mt). (52) 

For the pion-nucleon form factor we adopt the popular monopole parameterization 

G.^7v(t) = G^^iv~13, A^O.SGeV. (53) 

In the t — > limit {711^^ > 0) the Goldberger-Treiman relation shows up, which ties the axial-vector 
coupling to the pion-nucleon coupling gA = f-KG-^NN / \/^M [9l]. In the chiral limit (m,r — ^ 0) the 
axial-current is conserved and the r.h.s. of Eq. (152!) vanishes. Our ansaetze (149|) and (150|) are so 
designed that the resulting constraint 

Yim^^E^t) = H^,{t) with gp = AgAM^/GeY^ (54) 

is satisfied within the well-known relation of pseudoscalar and axial-vector coupling. If we now 
modify the relation between (partonic) pseudoscalar and axial-vector couplings by a small 
proportional term, 

^vai/sea ^ g^h^^ '^-^ + ^^' ^M^/GeV' , gA = gf + gT ^ 1-267 , (55) 
Z — a 



and set in the pion-nucleon form factor fl53|) the monopole mass A = ^J{1 — a)/a' , our model 
exactly satisfies the PCAC relation fl52|) . Moreover, within the generic Regge parameters fl36|) 
we find for the monopole mass the phenomenological value that is quoted in Eq. ( l53l) . Since we 
satisfy the PCAC relation and have a realistic description of the axial- vector form factor, see left 
panel in Fig. HI we also describe in the left panel of Fig. [5] the t-dependence of the pseudo-scalar 
form factor (IST!) measurements. 

We specify now the remaining parameter set {/?, 5/?, e} for both valence and sea quarks. As in 
Sect. 13. 2^ these parameters could be partly read off from large-x counting rules. Unfortunately, to 
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Figure 5: In the left panel the pseudo-scalar form factor is displayed: generic model (I50p prediction 
(solid curve), experimental [92] (triangles) and lattice (squares) data. In the right panel GPD model (jSOp 
predictions within parameters (j56p [thick] and ()57p [thin] and lattice results for generalized form factors 
are presented: j = 1 (dot-dashed, triangles) and j = 2 (dashed, squares). Same scale specification and 
lattice data as in Fig. [H 

our best knowledge, for E GPD such counting rules are not presented in the literature. We might 
expect that due to the non-conserved target helicity a E GPD decreases by two powers of (1 — x) 
faster at large x than a H GPD. Our E GPD Mellin moments fl50l) behave as for j ^ oo 

and, thus, in the limit x ^ 1 the zero-skewness E GPD vanishes as (1 — x)^"^^. For simplicity, we 
take this relative (1 — x) suppression factor in our model, i.e., we take for the E GPDs the same 
generic values for f3 (and also for 6/3) as for H GPDs, specified in Eq. (HSll . 

Now the e parameters can be used to tune the decrease of the E GPD Mellin moments (l50l) with 
growing j. To some extent they might be constrained by lattice data for generalized form factors 
(set VI of Ref. [SH]), shown in the right panel of Fig. O We use again our SU(6) symmetric valence 
scenario psl) within a flavor asymmetric sea and evolve our moments, specified at = 2 GeV^, to 
the scale //^ = 4 GeV^ of Ref. [89]. For j = 1 we find that e^'^' is slightly larger than one (triangles, 
thick dot-dashed curve). Since the j = 2 generalized form factor data (squares) are compatible 
with zero and very noisy, the sea quark parameter P*^*^ is essentially not constrained by data. We 
choose the parameter set 

gvalLat^^^^ gseag^n^^ ^^^^ 

which is fully compatible with lattice data (thick curves). Below we will also use the set 

gval = gsca _ Q ^ ^57^ 

which is, however, disfavored by the lattice measurement of the j = 1 generalized form factor, see 
thin dot-dashed curve. 
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Finally, we summarize our Ejj GPD model. It arises from plugging the PCAC relation fl55|) 
into the ansatz fl50l): 



^3)+ ^ ^3),val ^ ^3),sca ^ ^3)- ^ ^3),val ^ ^.^^ ^^g^ 

~(3) ^ AgAa' M\2-a- a^) (1 + j - a)^ B{1 + j - a,l + l3) - eBjl + j - a,l + f3 + 6f3) 
(2-a)(j-a^(t))(l+j-a(t))2 B{1 - a,l + (3) -eB{l ~ a,l + (3 + 6(3) 

the SU(6) symmetric valence scenario, cf. Eq. fl45l) . generic Regge trajectories (1361) . and by the e 
parameters (IS^ or, alternatively, (15 7p . 

3.4 The properties of a minimalist GPD model 

Let us shortly recall that at LO accuracy and for (nearly) fixed photon virtuality the GPD on 
the cross-over line rj = x entirely determines the 'CFFs' [9], see also Ref. [29] for a phenomeno- 
logical discussion to reveal this GPD from experimental measurements. From Eq. (1131) it follows 
that the imaginary part 

9m I ~ I (e, t, Q2) LO ^ 1 1 1 = ^ = g2) (59) 

of a 'CFF' is directly expressed by the GPD on the cross-over line. Moreover, for any GPD, i.e., 
this function has a double distribution representation, the unsubtracted 'dispersion relations' 

(e, t, Q') PV / dx -f^ { }{x,v = X, t, Q^) , (60) 




(^, t, Q2) LO ^ I (x, = X, t, Q^) , (61) 

offer an alternative way to evaluate the real part of CFFs from the GPDs on the cross-over line. To 
spell it out clearly, these 'dispersion relations' arise from the physical ones in the leading twist-two 
approximation, i.e., the spectral function and its support |,^| < .^cut/poi = 1 change [15], and they 
are in one-to-one correspondence with the LO approximation of CFFs f|T3|) . 

Only the skewness dependence of the GPD on the cross-over line is relevant for the LO de- 
scription of experimental data. To quantify this dependence, we define the r-ratio 

of the GPD on the cross-over line and on the line 77 = 0. If this ratio is one, we call it zero skewness 
effect. Note that the r-ratio is in a LO interpretation of DIS and hard exclusive electroproduction 
to some extent 'measurable' for H and H type GPDs. 
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To build a model for the full H^^^"^ and E^^^"^ GPDs, we relate them by a 'holographic principle' 
to their forward analog. In the following we prefer simplicity and take the minimalist 'holographic 
principle' in conformal Mellin space, see explanations in Ref. [49]. Namely, we equate the conformal 
GPD moments 

H,{r],t) = H,,{t)di'^,\r]) and 4(r/, t) = 4o(^) , (63) 

with the leading S0(3) PWs J = j + 1 and J = j, specified in Eqs. (l22i) and (123|) . respectively. 
The price of simplicity is that the GPD on the cross-over line is rigid. Our leading S0(3)-PW 
model fl63p can be equivalently set up in a minimalist "dual" parameterization or the Shuvaev 
transform [93] (where we understand that this name already implies the necessary restriction of 
the GPD transform [91]). The corresponding GPD and the double distribution have so far not 
been analytically evaluated, for other model examples see [111[95]. Numerically, the GPD can be 
evaluated within a Mellin Barnes integral [S]. 

We are here only interested on the phenomenological GPD aspect and need only the GPD 
on the cross-over line. Our specific model assumptions imply a characteristic small-x and large-x 
behavior, which can be read off from the imaginary part of the Mellin-Barnes integral (12T|) by 
means of Eq. (ISUj) . where we set = x. In the following we are not concerned about evolution 
and set = Ql in Eq. m\. 

The small-x behavior is determined by the leading Regge pole, which can be picked up by 
a shift of the integration path to the left hand side. One easily finds a partial factorization of 
t-dependence and a specific value for the r-ratio 



x-o ^ ' ' ^ r(3/2)r(2 + a(t)) ' x^o ^ ' ' ^ r(3/2)r(3 + a^(t)) ^ ' 

Since this is nothing but the Clebsch-Gordan coefficient in the conformal PW expansion fl2T]) . we 
call it for t = conformal ratic^ [19]. This ratio is only a function of the leading trajectory a{t) 
and it is valid for a restricted a{t) interval, i.e., —t is restricted by an upper bound. For a{t) > 
and a7r(t) > — 1 our models possess a positive skewness effect r > 1. 

The large-x behavior of a minimalist GPD on the cross-over line, can be derived from the large- 
j behavior of the Mellin-Barnes integrand (12T1) taken at large x. It arises from the Clebsch- Gordon 
coefficient, the skewness- zero GPD moments, and the behavior of SO (3) partial waves 

4y (^) - ^ r(5/2 + j)^^ ^'"''^^ + A = {0,1}. 



^The name is also motivated by the fact that a minimalist GPD model with leading Regge behavior in the 
small cc-region and for t = is equivalent to a conformal GPD model. Since effective Regge behavior is commonly 
present in a realistic GPD model, we name the ratio for shortness conformal. 
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Here we utilized the definitions (1221) and fl23|l in terms of liypergeometric functions and tlieir large 
j asymptotic for x ^ 1 [96]. Combining all together, and considering the Mellin-Barnes integral 
as a Fourier transform with conjugate variable ^/l — x"^ ~ V^V^ ~ x , we find that our minimalist 
GPDs behave at x — > 1 as 

H{x,x,t) ^ {l-x)^ (65) 
E{x,x,t) ^ (l-x)^ %^ ~(l-x)2. 

As required, our findings are consistent with the large-x result in the 'dual' modej^ [97]. The 
behavior (165|) also coincides with the large-x counting as it has been derived for the H and E 
GPDs from diagrammatic considerations [DP]- Hence, compared to the zero-skewness GPD, our 
GPDs at the cross-over line possess a large enhancement effect. We also note that in our models 



within linear Regge trajectories the t-dependence dies out at large x with y/1 — x. For fixed t the 
skewness ratio in the large x-region can be written as 

limr x,t 1 -a;"— ^ 7,/ , 66 

where F'^'^^ is a reduced GPD that encodes the t-dependence for both 77 = and t] = x cases. 
From the point of view of light-cone wave function modelling, |100l llOlj . double distribution 
models |lU2t I1U31 llU4j . and diagrammatical counting rules [HI] one expects that the t-dependence 
dies out faster for x — 1. It is here far beyond the scope to fix the presumable unrealistic feature 
in our models, e.g., by employing slightly non-linear Regge trajectories. 

For t = we illustrate the skewness effect in Fig. El where we show GPDs on the lines t] = x and 
?7 = as thick and thin curves, respectively. More specifically, we plot the dominant charge parity- 
even GPD models H^^^~^ (left) and xE^^^~^ (right), where the latter is rescaled by x. Thereby, the H 
pgj) and E fl58l) GPDs are specified by the parameter sets psj) and fl56l) . shown as dot-dot-dashed 
curve, and alternatively for the E GPD within the parameters flFTj) as dot-dashed curves. As one 
realizes in the right panel the E GPD models are qualitatively different, however, they approach 
each other in the small x-region. The characteristic small-x enhancement of a minimalist GPD is 
controlled by the conformal ratio (164]) . For our H and E GPDs it is given by r ^ 1.2 and r ^ 1.5, 
respectively. The large-x enhancement (1651) is clearly a big effect for the dot- dot- dashed models. 
We also observe that, compared to its forward analog, our minimalist GPDs on the cross-over line 
are enhanced in the valence region. Note that our E GPD on the cross-over line is always positive, 
while in the xQSM estimate it can be even negative [105] . 



^For H GPD one might utilize the S0(3)-PWs d^Q and so the power is (3/2 in agreement with Ref. [Wl [^ . 
Within S0(3)-PWs dp 1 polynomiahty is not completed and we obtain (/3 — l)/2 = 1 for /? = 3, which is consistent 
with the large a;-behavior in Radyushkin's double distribution ansatz with b ~ I [98j . 



26 



^ 3 



\ 




\ 




:\ \ 








\ \ 
\ \ 






■•-■■I 



3 



■ bj 2 



0.1 
X 















■ ^ \ 








1 





^. \ 


0.01 


0.1 


1 



Figure 6: The charge parity-even H'^^^'^ {x , x , t = 0) GPD (left) and scaled xE^^'^'^{x,x,t = 0) GPD 
(right) with parameter sets (|45|) and (|57|) are shown as dot-dot-dashed curves and for parameters (|56|) as 
dot-dashed. The corresponding GPDs at = are displayed as thin curves. 

Let us now point out the particularities of the charge parity-even GPD. This GPD steeply 
rises as with decreasing x, see the flat behavior of the scaled xE GPD in the right panel 

of Fig. O Since the small intercept = —fnla' is negative, it is guaranteed that for t < the 
lowest Mellin moment of this GPD exists and that an unsubtracted dispersion relation fl^Ul) is 
sufficient. Because of this singular behavior at x = 0, the charge parity-even 'CFF' (!60l) possesses 
a large real part that can be loosely viewed as a 1/^ proportional constant. To isolate this term, 
we introduce an over subtraction in the dispersion relation (1601) at the point ^ = oo: 



3f?e£(3)+(e,t) 



1, 



X- 



-E(=^)+(x,x,t) +iQ(t), Ceit) 



1 _ 

2 / dx£^^'^+{x,x,t) . (67) 





'0 e 

The subtraction constant in this new 'dispersion relation' is evaluated from the GPD on the cross- 
over line and, thus, it encodes the skewness effect. In the vicinity of t = m^, we find by means of 
the Mellin-Barnes integral ( l2T|) the pion pole: 

Lo 3 gp 



Cf-(t) 



non pion-pole contributions . 



2ml-t 

This is nothing but the result in the parameterization of Ref. p3l [T6 
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c(^' - |C'^poic(^) 5 C'^-poic(^) - J dx ^_ ^Ji _ 
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where the residue is provided by the asymptotic pion DA 

/^6u(l - u). 



U 



(69) 



(70) 



Few comments are in order, which are partially adopted from discussions of parity-even GPDs 
and the so-called D-term [IDH] in Refs. ^ [29]: 
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i. In the parameterization of Ref. [HI [16] the pion-pole hves only in the central GPD region and 
it is decoupled from the smooth GPD part. This might be the natural answer that arises from a 
dynamical model calculation, in particular, if the chiral limit is taken. In such a model the cross 
talk of imaginary and real part is interrupted. On the other hand we can consider the chiral limit 
of the subtracted dispersion relation (1671) . Although the intercept is now = —1, the integral 
still exists for t < and for t > the answer can be obtained by an analytic regularization of the 
integral. 

a. One might also reverse the logic, and utilize the subtraction constant (1671) as a sum rule that 
constraints the GPD on the cross-over line. 

in. The residue of the pion pole in the minimalist GPD model, see Eq. fl68p . is equal to that one 
in Eq. fl69|) . obtained from the asymptotic DA. This is not accidental and it can be viewed from 
a more general perspective, reflecting the crossing relation between t- and s-channel. Namely, 
a family of asymptotic generalized DAs, labelled by the t-channel angular momentum J is the 
crossing analog of the minimalist GPD. 
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Figure 7: Left panel: ?R.e£^ {£,,t) /£^_p^^^{(^,t) ratio, cf. Eq. ([69]) . versus ^ for fixed t = (thin) and 
t = — 0.2GeV^ (thick). Right panel: cross section ratio daL{xB,t'\H , E) /daL{xB,t'\H = 0, E) versus xb 
for t' = (thin) and t' = -0.2 GeV^ (thick) with = 2GeY'^. The dot-dot-dashed and dot-dashed 
curves refer to the parameterizations (j56|) and ()57p . respectively. 



In Fig. [7] we illustrate the role of the pion pole contribution for GPD modelling. In the left 
panel we show the ratio of the real part of our 'OFF' model to that one of Ref. [15], entirely given 
by the pion pole contribution fl69|) . versus ^. For t = the ratio is shown as thin curves and as 
expected it is about one. It slightly gets smaller with decreasing ^, where the comparison of the 
dot-dot-dashed fISBl) and dot-dashed curves reveals that the parameterization of the GPD on 
the cross-over line has some influence. For non-vanishing —t the difference is more pronounced. 
Since we neglected the residual t-dependence in the pion pole ( 169!) . the ratio considerably decreases 
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for -t = 0.2 GeV^ (thick curves). In the right panel we display the ratio of the unpolarized cross 
sections ([3]) for the full H and E GPD model fl49|58p to that one in which H is neglected. Thereby, 
we set t = tyain + t' . Clearly, only for t' = and in some limited xb interval this ratio is around one 
and so the E dominance is justified. We conclude that in general kinematics, even in the valence 
quark region for smaller values of —t', the hard exclusive 7r~^ electroproduction cross section might 
not be dominated by the pion pole (!69l) . 



4 Measurements versus formalism and minimalist model 

We like now to confront our minimalist GPD models with the exclusive vr"*" electroproduction cross 
section [30] and single transverse proton spin asymmetry [36] measurements from the HERMES 
collaboration and the longitudinal photoproduction cross section, measured in Hall C at JLAB, [211 
[32] . However, first we will define the observables that can be related to the coUinear factorization 
approach. We like then to comment on the theoretical problems of the collinear factorization 
approach and we suggest a simple prescription to cure its LO approximation. We finally rely on 
this prescription and compare our ad hoc GPD models and a tuned one with experimental data. 

It is common to convert the electroproduction cross section to the virtual photoproduction one 
by means of the Hand convention for the photon fiux factor Fy: 

da^'P-"^^ 1 /-^^ da-^ ^ ttem 1-XB 

"■^37=^73 — 3777' = -— (71) 



dt Tv{xb,Q)Jo dQ^dxBdtd<p' 2ti I - e xbQ^ 

The photoproduction cross section can be separated in the transverse and longitudinal part 

_ daj:'^^""^ ^ ^ daf ^^""^ ^^^^ 
dt dt dt 

Utilizing the Hand convention fl7T]) . we can read off the perturbative prediction for the longitudinal 
polarized photoproduction cross section from Eq. ([2]): 

'"^^""^'^-^ /a-aiH..r-74ieE..p-2e^^eH..E:, 



dt ^/iT^Q'^l-XBy - ^' - I 4M2' 



(73) 



If the extraction of the longitudinal photoproduction cross section is experimentally not reach- 
able, a comparison of GPD model predictions for the longitudinal cross section fl73|) and elec- 
troproduction cross section measurements requires a quantitative understanding of its transverse 
part, too. Utilizing the R ratio, we might express the measured differential photoproduction cross 
section (fTTj) to its longitudinal part (fTSj) as follows: 



dt 



1 + 



s(a;B,Q2)i?(xB,t, Q2) 



e(xB,Q^)^^S ' R = -^^ 7- (74) 
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Note that the factorization theorem [T] tells us only that the inverse ratio 1/ R formally vanishes 
with 1/ Q^, however, it does not state that it is numerically small at accessible photon virtualities. 
Rather it should be expected that l/i? is numerically enhanced by non-factorizable contributions, 
which might also include rather large logarithmical modification of its expected 1/ fall-off; for 
an example see the perturbative study to the proton Pauli form factor in Ref. [107j . Indeed, 
measurements in Hall C at JLAB reveal that l/i? is rather sizeable for < 4 GeV^ and W ~ 
2GeV [Ml [32]. 

Unfortunately, to our best knowledge, also a model understanding of the 1/R ratio is at 
present not fully reached. For instance, the Regge-inspired model [33] describes at such low W 
the longitudinal and the real transverse photoproduction cross section measurements, however, 
underestimates the virtual transverse one. On the other hand it is in fair agreement with the 
HERMES measurements for the virtual photoproduction cross section ^U]. A refined version of 
this model, tuned to low energy data from JLAB [31], does describe most of HERMES data, too, 
however not those at large photon virtualities and relative low energy [35] . In the models the pion 
pole is reggeized and it essentially contributes only to the longitudinal cross section, which for 
\t'\ < 0.5 GeV^ is dominant. Hence, for smaller —t' values both models suggest that we can set 
to zero. 

Only the sin{(f) — (ps) harmonic of the single transverse proton spin asymmetry 

— — = Aut"^ sm{(j) — (j)s) + other harmonics , (75) 

is counted as a further twist-two observable [16], while other harmonics are at least formally 
suppressed by 1/Q. In the notation of Ref. [108] the asymmetry is written in terms of so-called 
polarized photoabsorption cross sections a. 



V ■ 

mn' 
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e cr, 



00 



where the subscripts (superscripts) refer to the photon (target) helicity. Both terms and 
aijl^ + oc daT arise from the exchange of a transverse polarized photon and, thus, are counted 
as power suppressed non-factorizable contributions. We neglect both of them and from the 
longitudinal cross section ([3]) we read off the twist-two approximation for the asymmetry 



^sin{<l>-^s) Tw-2 \/ ^ TT 



(1 - e)|H.+ P - ji^l^E^A' - 2e5JeH,+E;, 



(77) 



Note that the authors of Ref. [16] argued that non-perturbative contaminations mostly cancel in 
this asymmetry and that a cancellation of perturbative corrections has been numerically estab- 
lished within some GPD models in Refs. [T8l [25] . 
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As discussed in Sect. [H it is an intricate task to judge on the coUinear factorization approach, 
apphed to the hard exclusive meson electroproduction. We recall that the hard amplitude is rather 
analogous as for the pion form factor, where the momentum fractions of the incoming minimal pion 
Fock state are replaced by those of the emitted and absorbed quarks [TS], see Fig. [T]d. Moreover, 
the pion pole plays an important role at very small —t and not too small x-q. Thus, the theoretical 
descriptions of the hard exclusive vr"*" electroproduction and pion form factor have similarities. 
The perturbative framework for the pion form factor has been widely discussed in the literature, 
for critics see Ref. |1U91 IllU] . and various improvements were suggested. One might use these 
recipes to improve the description of hard exclusive meson electroproduction, too, e.g., going 
back to the overlap representation of Drell and Yan in which the meson wave functions contains 
transverse degrees of freedom [88], modifying the factorization by incorporating the Sudakov 
suppression [lll^ I112j . or a scale setting prescription that reshuffles perturbative corrections to 
the non-perturbative region of the coupling [113] . 

A popular GPD code is called VGG plus "power corrections". Here the "power corrections" 
result from a 'hand-bag' calculation in which transverse degrees of freedom (and only the pion 
pole in the E GPD) are taken into account [13]. It should be clear from the very beginning 
that these "power corrections" can not be considered as autonomous contributions and that their 
interplay with higher Fock state components or so-called soft contributions is not obvious; see 
also for illustration the light-cone sum rule analysis in Ref. |114] . In some sense this recipe is 
analog to the Drell- Yan formula for the pion form factor and here the conceptual problem of 
such an approach is obvious. Namely, the Drell- Yan formula is the starting point and integrating 
out the transverse degrees of freedom yields the coUinear factorization approach. This procedure 
generates perturbative corrections in the leading twist-two sector [5] and also power suppressed 
contributions. Strictly spoken, even perturbative and power suppressed contributions can not be 
autonomously considered as they have a cross talk. Already, for some time it was observed that the 
sign and magnitude of "power corrections" [15] and perturbative NLO corrections in both DVCS 
|1151 H5] and hard exclusive meson electroproduction [18] match. Although this observation is not 
quantitatively understood, it reveals the problem of double counting. Hence, one can not take the 
notion VGG plus "power corrections" literally. 

It would be of immense phenomenological interest to estimate the size of perturbative and 
non-perturbative corrections to the LO formulae ([6]), ([8]) of the coUinear factorization approach. 
Since not much is known about power suppressed corrections, we will here assume that they can 
be neglected. The biggest portion of the large perturbative NLO corrections is proportional to the 
first coefficient of the (3{as) function [18], which governs the scale change of the running coupling 
([9]). These large corrections can be viewed as an indication for a break down of perturbation 
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theory. This calls for a resummation procedure, which naturally induces its own ambiguity. Since 
the LO approximation of the hard scattering amplitude ([8]) is proportional to the running coupling, 
a big portion of the NLO corrections can be absorbed by changing the renormalization scale fiR. 
The ambiguity we fix by the Brodsky-Lepage-Mackenzie (BLM) scale setting prescription |116j 
(for the treatment of complex amplitudes see Ref. |117j ). yielding 

a.(/z^L^ = exp{-/(0/2}Q) with e-^(«)<0.1. 

The resulting renormalization scale, i.e., the function f{^), depends on the GPD itself, however, 
it will be certainly driven into the infrared region and so the coupling itself is now plugged by 
non-perturbative physics. 

The BLM scale in hard exclusive electroproduction is rather small [TSl I118j . similar as in 
the perturbative description of the pion form factor |113l I119j . In the following we do not aim 
for an "optimal" description of the hard exclusive 7r+ production rather for a semi-quantitative 
understanding. Therefore, we simply assume that the effective couplings, arising from the BLM 
scale setting procedure and the freezing in the infrared region, have in both processes the same 
value. Hence, within a very low renormalization scale we can employ the LO approximation for 
the (scaled) pion form factor, 

Q'FAQ')^^87raff^ll{Q'), (78) 

to find the value of this effective coupling. As mentioned, the inverse moment X^r, defined in 
Eq. fllOj) . can be directly revealed at LO from the photon-to-pion transition form factor. Its value, 
determined from CLEO data [120j . would be slightly smaller than one, see, e.g., Refs. [1211 1122"] 
for comprehensive discussions. We set it here generically to one, i.e., X7r(Q^) = 1. This value also 
arises from the asymptotic pion DA fl70|) however, also pion DAs with a more complex shape might 
approximately possess the same value [123^ 11241 1122j . For the scaled pion form factor (I78p we 
take the value 0.36, which is for > 1.5 GeV^ and within 1-cr standard deviation compatible 
with the model dependent extraction of the form factor from the vr"*" electroproduction |125j . and 
so we obtain the following value for the effective coupling: 

af = 0.8 . (79) 

Both the photon-to-pion transition and the elastic pion form factors are described at LO. Note 
that also within the BLM setting prescription sizeable NLO corrections, related to the Sudakov 
suppression, affect the LO description of the pion form factor |113j . 

To provide predictions from our minimalist GPD models, we rely now on the LO approxima- 
tion of the TFFs 0151) . where the effective coupling is specified in Eq. fl7^ . Thereby, we employ the 
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Figure 8: HERMES measurements [30] of the differential cross section da'^*^^'^'^^ /dt' versus minimalist 
GPD model predictions: only with pion pole contribution (dotted), including H (dashed), full model 
with the setting (j57p and frozen coupling (dot-dashed), same as dot-dashed within parameters (j56p (dot- 
dot-dashed), and versus a tuned GPD model ([80]) (solid). 

Mellin-Barnes integral (1211) and the forward Mellin moments ( l63l) . which are skewed by the PW 
amplitudes (149|) and (l58l) . We like to add that the integration contour must be chosen carefully 
and that the dispersion relations (160|) and (16T1) served us to check our numerics. In the following we 
list the various predictions for the differential cross section ([71]) that are displayed together with 
the HERMES measurements [3D] in Fig. [S] where statistical and systematical errors are added in 
quadrature (same applies for Fig. [9] and Fig. [TO]) . In the confrontation of model predictions and 
measurements we should bear in mind that at small —t' < 0.5 GeV^ the longitudinal cross section 
should dominate and that the large —t' region, in which the neglected transverse cross section 
plays a role, is a priory not suited for the coUinear factorization approach. 



dotted: pion pole and H = 0, running coupling 

We recall first the model prediction that relies on an overwhelming role of the pion pole, shown 
in Fig. 4 of Ref. [SU] as dot-dashed curve. We imitate this model by setting H GPD to zero 
and replacing S by its real part within the parameter set (1571) . see also left panel in Fig. [3 The 
resulting cross section is plotted in Fig. [8] as dotted line, which resembles the corresponding curve 
in Fig. 4 of Ref. [30]. As one realizes the model fails to describe the normalization of the cross 
section at very low —{t') ~ 0.08 GeV^, where the discrepancy can reach one order of magnitude 
at larger values. Since Regge-inspired model estimates [331 [3S|, e.g., compare dashed and 
dotted curves in Fig. 4 of Ref. [30], state that the 1/R ratio should be close to zero at t' = 0, a 
popular conclusion drawn from the observation is that so-called "power contributions" are needed. 

dashed: pion pole and minimalist H model, running coupling 
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We take now the same model for £ as above and include the H GPD within the parameter set 



[5|) . displayed as dashed line. As expected from our considerations in Sect. 13. 4[ the inclusion of 
H yields a flatter t-dependence, compare dashed with dotted curves. One also realizes that for the 
lowest xb values (left panel), which means also low = l.4GeV^ the cross section at small —t' 
increase by roughly a factor of two and so the discrepancy of model estimate and measurement 
diminish, see also the left panel in Fig. [3 However, the discrepancy increases with growing 
and reaches approximately one order of magnitude at (Q^) = 5.4 GeV^. Since the LO prediction 
of the cross section is proportional to the square of the running coupling ([H]), its normalization 
changes within the canonical scale setting prescription for the renormalization scale /i/? = Q by 

^^7*P-n-+(a;g^t, = 5.4 GeV^) al{Q^ = 5.4 GeV^) _ ^ ^ 
da^P^^^+ixB^t, Q2 = 1.4 GeV^) ^ a^^{Q^ = 1.4 GeV^) 

Hence, we might be tempted to conclude that half of the discrepancy, observed at = 5.4 GeV^, 
arises from logarithmic scaling violation, mainly induced by the running coupling. 



dot-dashed: pion trajectory E model ( [57| j and minimalist H model, frozen coupling 

If we now freeze the coupling (1791) . our model with reggeized pion pole (SJrnT 7^ 0), describes data 
rather well at smaller —t' values, see dot-dashed curve. The growing discrepancy within increasing 
— t' might be understood as an increase of the transverse cross section. Indeed, the dot-dashed 
curve is rather close to the Regge model prediction |33] for the longitudinal cross section, shown in 
Fig. 4 of Ref. [30] as dotted curve, while the prediction for the full differential cross section, shown 
there as dashed curve, matches the experimental measurements. Hence, if one would borrow the 
1/R ratio from the Regge-inspired model, we would presumably conclude that our GPD model 
describes data over an astonishing large \t'\ range. 

dot-dot-dashed: pion trajectory E model ^5R) and minimalist H, frozen coupling 

Unfortunately, the parameter set ( l57j) is rather disfavored by lattice data, see dashed curve in the 
right panel of Fig. [5l Moreover, as we will see below in Fig. [101 it predicts that the single trans- 
verse proton spin asymmetry is of order —1 and so it might be considered as ruled out. Taking 
the parameter set fl56l) . which is more favored by lattice simulations, we obtain the dot-dot-dashed 
curve in Fig. [HI Now the discrepancy of the normalization at small t' increases again with growing 
Note also that the t'-dependence is now flatter and if we would correct the normalization for 
large value of in the right panel, the model would perfectly describe data. In other words our 
model description would then state that the cross section is dominated by the longitudinal part 
of the cross section over a large lever arm in —t. This would then contradict the Regge-inspired 
model predictions 
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Figure 9: Total cross section fj7*p^"T^ versus measurements from HERMES [3D] are sliown together 
with a tuned (solid) and various ad hoc GPD model estimates with 1/R = 0. Same models as in Fig. [8l 



solid curve: pion trajectory E model (51) and tuning H, frozen coupling 

As we demonstrated, our minimalist GPD models provide us with a variety of predictions and 
confronting them with experimental data might imply contradictory conclusions about the failure 
or successes of the collinear factorization approach. Hence, as advocated in Ref. [29], it is more ap- 
propriate to change the GPD approach. We now assume that the collinear factorization approach 
with a frozen coupling works and ask now for the functional form of the GPD at its cross-over 
trajectory, which entirely determines the imaginary and real part of CFFs to LO accuracy. Instead 
of introducing a new parameterization for GPDs on the cross-over line we employ here the Mellin- 
Barnes integral again. Only for illustration we consider it as representation for a complex valued 
function, satisfying the dispersion relations fl60|) and fl6T]) . where the physical/partonic meaning of 
the parameters is mostly lost. From an "eye fit", where we already included preliminary HERMES 
data for the single transverse proton spin, we find that the parameters set 



^vai "fit" Q_gg ^ ^vai "fit" , 5/5^^1 "=" 1/2 , /i^^^ "=" 1 , "=" 5 , 6(3''''' "=" 1 , (80) 
F^i , pf =" 1 , =" , "=" 3 

will lead to a good description of data. In Fig. [S] we realize that the solid curve indeed describes 
the small —t' data points and leaves some space for the contribution of the transverse cross section. 

In Fig. [9] we also display our model results for the total cross section cr7*p^"'r+ versus Q^. Such 
confrontation is sometimes viewed as a test for the onset of the perturbative approach. However, 
we should also spell out some warnings. First at all, the total cross section o"7*p^a^a/ ^ 
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Figure 10: Preliminary HERMES measurements [36] of the single transverse proton spin asymmetry 
^sin(</) <j>s) ygj,g^g tuned (solid) and various ad hoc GPD model predictions, same as in Fig. El 



pure twist-two quantity and secondly even if we would have a measurement of its longitudinal 
part cr^ ^^^^ it contains necessarily contributions that arise from phase space regions for which 
the coUinear factorization approach is perhaps not applicable. Moreover, since the phase space, 
accessible in experiment, is restricted, the kinematical variables might be even correlated. Naively, 
one would expect a scaling and if it would show up in the data one might say that the onset 
of scaling justifies the perturbative approach. Also the reversed logic has been often stressed in 
the literature. Certainly, in the first place the running coupling, appearing in the hard amplitude, 
and also logarithmic scaling violation, predicted by evolution, demolishes this logic and reminds us 
that logarithmic scaling violations are part of the perturbative framework. Indeed, our two models 
within running coupling (dotted and dashed) have a steeper Q^-dependence as the models that 
were treated with a frozen coupling (dot-dashed, dot- dot- dashed, and solid) and dashed). The 
slope of our "fitted" GPD model prediction (solid) fairly follows the data. That it underestimates 
the total cross section measurement is natural, since we describe only the differential cross section 
at smaller —t' values. The differences might arise from the transverse cross section. 

We now confront our model predictions in Fig. [10] with a preliminary measurement of the 
single transverse proton spin asymmetry Aux*^ which we describe within the approximation 



7]). We realize from this formula that the asymmetry is proportional to 
j^H^-^s) Tw-2 _ (^^^ud\^ (^cjmH"'^) + (Sm^"-^) (^eTC" 

For the pure pion pole model (dotted) the asymmetry trivially vanishes. In a pion pole model 
with non-vanishing H GPD the asymmetry, shown as dashed curve, is usually predicted to be 
negative with a large modulus that is rather fiat with respect to the t'-dependence; except that it 
due to the kinematical factor A/t'/^mm in the asymmetry (1771) steeply vanishes at the kinematical 
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boundary. These characteristic features arise from the fact that the imaginary part of the Ti. CFF, 
equivalently, the H GPD on the cross-over hne, is in all popular GPD models positive and the 
imaginary part of S plays no or not an essential role, see Eq. ( ISTI) . Obviously, in a twist-two 
interpretation such standard GPD models [121 [IB] are ruled out from the preliminary HERMES 
measurement |36j . 

The dot-dot-dashed curve illustrate that the prediction from the reggeized pion pole GPD 
model in the version (156!) looses the characteristic t'-dependence and it can even possess a sign 
change. In this GPD model the real part of the CFF S becomes at larger values of —t' in the 
vicinity of xb ~ 0.1 negative, while Q'mTi remains positive. However, also this twist-two model 
prediction is still disfavored by experimental data. 

If we now give up our minimalist model for the GPD H and employ the pragmatic parame- 
terization fISU]) . the resulting solid curves are in fair agreement with the experimental points. In 
this model the real part of S is also for larger —t' values positive (in contrast to the dot-dot- 
dashed model, since we have now e^^^ = e^^^ = 0), while now the imaginary part of Ti. changes at 
-t' ~ 0.3 GeV^ the sign and so also the asymmetry. Note also that both terms on the r.h.s. of 
Eq. flSip are sizeable, where the second term behaves differently than the first one: it is positive 
at small —t' values and becomes negative at larger —t' values. Thus, also delicate cancellations 
appear in the asymmetry prediction, which finally determine the shape of the solid curve. 

From the confrontation of our GPD models with the transverse target spin asymmetry mea- 
surements we conclude that a "predictive power GPD model approach", as it has been emphasized 
for the overwhelming role of the pion pole, might be misleading. The phase difference between the 
'CFF' Ti and S are not predictable at larger values of —t'. As said before, relying on the twist-two 
dominance of this observable one can "fit" experimental data to reveal GPDs. Unfortunately, we 
have also illustrated that rather different GPD model scenarios are potentially able to describe 
the HERMES measurements. 

Finally, we confront in Fig.[TT]our models with the measurements of the longitudinal differential 
cross section daj^ ^^^'^'^ /dt from Hall C collaboration [32], where we do not account for systematical 
errors. As one realizes also this low energy data are described by the GPD model fl80l) . "fitted" to 
HERMES data, while all other models underestimate the measurement. In particular, the pion 
pole model (dotted) with running coupling fails in the normalization and it possesses a too flat 
t-dependence. Freezing the coupling, the cross section estimate will be roughly enhanced by a 
factor five and the pion pole model without H is becoming closer to experimental data. This 
estimate is similar to our dash-dotted curve, based on the E parameterization fl57j) . However, we 
see from the alternative E parameterization fl56|) . shown as dot-dot-dashed curve, that the model 
predictions are sensitive to the shape of E GPD, too. Now the model roughly underestimates 
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Figure 11: Measurements of the longitudinal differential cross section daj^ p^i^'^ Hall C at JLAB 

[32] versus various GPD model predictions, same as in Fig. [8l Only statistical errors are included. 

experimental data by a factor three. Employing the parameter set fISD]) . we enhance the GPDs 
in the large x region and describe experimental data. This enhancement feature of cross sections 
at low energy has been also observed in the hard p electroproduction and it was modelled within 
a D-like addenda [26]. We like to emphasize the intrinsic duality aspect of GPDs that such an 
enhancement is naturally tied to the GPD behavior on the cross-over line, see also the spectator 
model [104] . We add that our "fitted" GPD model also describes HALL C data from Ref. [3T] 
and that the VGG plus "power corrections" code [T5| is also in fair agreement with HALL C data, 
too. 

5 Summary and conclusions 

In this paper we have applied generic GPD modelling |5U] to set up zero-skewness H and E GPDs 
in the iso-vector sector. Thereby, we were confronted with the problem that Regge phenomenology 
within unnatural parity exchanges is challenging. Relying on the restoration of chiral symmetry 
in hard or even high-energy processes, we conjectured a master trajectory for odd J"*""*" exchanges. 
Such a restoration might be also somehow indicated by the appearance of the 03 (1875) meson, 
which is listed by the PDG under further states. The hypothetic trajectory naturally explains 
the rise of the polarized DIS structure function gi at smaller x-q- We illustratively took a SU(6) 
symmetric valence scenario for the net contribution. In such a scenario the 'spin puzzle' in the iso- 
vector valence quark sector is trivially resolved. The incorporation of t-dependence in our generic 
GPD model is based on a Regge pole product ansatz in Mellin space, where the PCAC Ward- 
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identity could be satisfied within a linear pion trajectory. Both form factor and generalized form 
factor lattice measurements are then well described. Certainly, the model can be refined, were in 
accordance with spectroscopic data a non-linear pion trajectory should be taken into account. 

Generic arguments to skew a GPD and to find so its value on the cross-over line are only 
given by large-x counting rules. Compared to the skewness-zero GPD, they are enhanced in the 
large-x region. Moreover, the skewness dependence of GPDs could be also constrained by precise 
lattice measurements, expected to be available somewhen in future, however, it is clear that this 
information is not sufficient to pin down a GPD on the cross-over line. Here we set up ad hoc 
models within the minimalist skewing prescription in which only the leading SO (3) partial wave 
in the expansion of conformal moments is taken into account. Our GPD models on the cross-over 
line possess the expected large-a; behavior, however, the suppression of t-dependence for x ^ 1 
might be too weak. More flexible models can be built by including non-leading SO (3) partial 
waves and they can be tuned to experimental data |49j . 

To confront the ad hoc GPD models with experimental measurements, we relied on the collinear 
factorization approach to LO accuracy. Thereby, perturbative corrections were taken qualitatively 
into account by incorporating them in the coupling via the BLM scale setting prescription. This 
might to some extent cure the perturbative approach. However, the coupling itself is now consid- 
ered as a non-perturbative expansion parameter that must be fixed from experimental data, too. 
In such a procedure it is rather likely that the coupling looses its universality and it should be 
rather considered as an effective one. Nevertheless, we end up with a semi-quantitative framework 
for the GPD phenomenology that can be tested versus experimental measurements. 

We demonstrated then that the hard tt"*^ production cross section for HERMES kinematics at 
small —t' can be described within our minimalist GPD models. Thereby, the H GPD is important, 
too, and so the pion pole in the E GPD plays an essential, however, not an overwhelming role. 
We also considered the single transverse proton spin asymmetry, preliminarily measured by the 
HERMES collaboration, where we assumed that power suppressed contaminations, induced by 
the exchange of transverse polarized photons, are small. For small —t' values the sign of the 
asymmetry is determined by both the pion pole and the sign of the H GPD on the cross-over line. 
In our ad hoc models the sign of both contributions is positive, in agreement with the experimental 
findings. However, or ad-hoc models are now disfavored for larger —t' values, where they predict 
as other popular models a large negative asymmetry. Tuning the GPDs on the cross-over line, we 
were able to describe the measurements, too. For larger —t' the pion pole is not overwhelming 
and the relative phase between 7i and S 'CFFs' has an intricate variable dependence, which can 
hardly be predicted. In our perturbative framework the same ad hoc models are also disfavored 
by the longitudinal photoproduction cross section measurements in HALL C at JLAB. 
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Our phenomenological studies illuminated the problems in the extraction of GPDs from hard 
exclusive meson electroproduction measurements. Finally, let us clearly spell out the challenges. 
On one hand GPDs are intricate functions, depending on the momentum fraction, the skewness 
parameter, and the momentum transfer squared. Thus, one can not expect that popular GPD 
models, which we mainly consider to be ad hoc, are able to describe experimental data. This 
proble nE can be resolved within flexible GPD models used in robust numerical fitting routines 
[inillQ]. So far this approach has been applied to deeply virtual Compton scattering in the small- 
xb region and within 'dispersion relation' technique for fixed target kinematics [H]. Indeed, such 
an approach allowed for the first time to access the H GPD at leading order and revealed so the 
failure of ad hoc models, used before. We emphasize that for the hard exclusive meson electropro- 
duction the theoretical framework is not well understood. To judge on the coUinear factorization 
approach and its improvements, one can utilize the universality of GPDs. Namely, from a widely 
phenomenological application of the GPD approach within flexible GPD models or 'dispersion 
relations' one might test its internal consistency. If such a phenomenological consistency check is 
successful, one simultaneously reveals the GPDs, mainly on the cross-over line. 
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First note added 

During the preparation of the manuscript we were becoming aware of (preliminary) photon-to- 
pion transition form factor measurements from the BaBar collaboration, which are now available 
in Ref. [128j . These experimental data do not support the canonical scaling, which we have im- 
plicitly assumed, to fix the inverse moment of the pion distribution amplitude. They challenge our 
interpretation of the connection among photon-to-pion form factor, pion form factor, and hard 
exclusive vr"*" electroproduction, however, the normalization of our numerical results is essentially 
adopted from the pion form factor. Nevertheless, the BaBar measurement [128j yield a 'crisis' 
for the perturbative description of pion related processes and the present understanding of the 

^Alternatively, one might first ask for the value of CFFs by utilizing analytic formulae [TUl SI] or fitter codes 
|1261 1127] . where the later method is accomplished with error estimates. 
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partonic substructure of the pion. First comments and interpretations of the BaBar measurement 
are given in Refs. [122[ I129[ 1130] . Since the pion-to-photon transition form factor is a key process 
for a QCD understanding, we would hke to encourage the BELLE collaboration to cross check 
the measurement and we recall that a double tagged experiment might provide a deeper insight 
into the 'pion crisis'. We add that the theoretical description of all exclusive B decay into pion 
channels is affected by this 'crisis'. 

Second note added 

After our manuscript was essentially finished, S.V. Goloskokov and P. KroU submitted their work 
An attempt to understand exclusive vr"*" electroproduction |131j . The authors rely on a 'hand-bag' 
approach, which can be considered as a model improved coUinear GPD approach. The advantage 
of such an approach is that various observables can be successfully described, which provides 
some trust in the underlying GPD framework. We try to stay as close as possible to the coUinear 
factorization approach in the hope that GPD universality and a systematical improvement can 
be ensured in a maximal manner. This we consider as a requirement that allows aiming for the 
primary goal, namely, to reveal GPDs from experimental measurements. 
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